Water buffalo genome science comes of age.

The water buffalo is vital to the lives of small farmers and to the economy of many countries worldwide. Not only are they draught animals, but they are also a source of meat, horns, skin and particularly the rich and precious milk that may be converted to creams, butter, yogurt and many cheeses. Genome analysis of water buffalo has advanced significantly in recent years. This review focuses on currently available genome resources in water buffalo in terms of cytogenetic characterization, whole genome mapping and next generation sequencing. No doubt, these resources indicate that genome science comes of age in the species and will provide knowledge and technologies to help optimize production potential, reproduction efficiency, product quality, nutritional value and resistance to diseases. As water buffalo and domestic cattle, both members of the Bovidae family, are closely related, the vast amount of cattle genetic/genomic resources might serve as shortcuts for the buffalo community to further advance genome science and biotechnologies in the species.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle.

BACKGROUND: The recently constructed river buffalo whole-genome radiation hybrid panel (BBURH5000) has already been used to generate preliminary radiation hybrid (RH) maps for several chromosomes, and buffalo-bovine comparative chromosome maps have been constructed. Here, we present the first-generation whole genome RH map (WG-RH) of the river buffalo generated from cattle-derived markers. The RH maps aligned to bovine genome sequence assembly Btau_4.0, providing valuable comparative mapping information for both species. RESULTS: A total of 3990 markers were typed on the BBURH5000 panel, of which 3072 were cattle derived SNPs. The remaining 918 were classified as cattle sequence tagged site (STS), including coding genes, ESTs, and microsatellites. Average retention frequency per chromosome was 27.3% calculated with 3093 scorable markers distributed in 43 linkage groups covering all autosomes (24) and the X chromosomes at a LOD >or= 8. The estimated total length of the WG-RH map is 36,933 cR5000. Fewer than 15% of the markers (472) could not be placed within any linkage group at a LOD score >or= 8. Linkage group order for each chromosome was determined by incorporation of markers previously assigned by FISH and by alignment with the bovine genome sequence assembly (Btau_4.0). CONCLUSION: We obtained radiation hybrid chromosome maps for the entire river buffalo genome based on cattle-derived markers. The alignments of our RH maps to the current bovine genome sequence assembly (Btau_4.0) indicate regions of possible rearrangements between the chromosomes of both species. The river buffalo represents an important agricultural species whose genetic improvement has lagged behind other species due to limited prior genomic characterization. We present the first-generation RH map which provides a more extensive resource for positional candidate cloning of genes associated with complex traits and also for large-scale physical mapping of the river buffalo genome.

Application of dissociation curve analysis to radiation hybrid panel marker scoring: generation of a map of river buffalo (B. bubalis) chromosome 20.

BACKGROUND: Fluorescence of dyes bound to double-stranded PCR products has been utilized extensively in various real-time quantitative PCR applications, including post-amplification dissociation curve analysis, or differentiation of amplicon length or sequence composition. Despite the current era of whole-genome sequencing, mapping tools such as radiation hybrid DNA panels remain useful aids for sequence assembly, focused resequencing efforts, and for building physical maps of species that have not yet been sequenced. For placement of specific, individual genes or markers on a map, low-throughput methods remain commonplace. Typically, PCR amplification of DNA from each panel cell line is followed by gel electrophoresis and scoring of each clone for the presence or absence of PCR product. To improve sensitivity and efficiency of radiation hybrid panel analysis in comparison to gel-based methods, we adapted fluorescence-based real-time PCR and dissociation curve analysis for use as a novel scoring method. RESULTS: As proof of principle for this dissociation curve method, we generated new maps of river buffalo (Bubalus bubalis) chromosome 20 by both dissociation curve analysis and conventional marker scoring. We also obtained sequence data to augment dissociation curve results. Few genes have been previously mapped to buffalo chromosome 20, and sequence detail is limited, so 65 markers were screened from the orthologous chromosome of domestic cattle. Thirty bovine markers (46%) were suitable as cross-species markers for dissociation curve analysis in the buffalo radiation hybrid panel under a standard protocol, compared to 25 markers suitable for conventional typing. Computational analysis placed 27 markers on a chromosome map generated by the new method, while the gel-based approach produced only 20 mapped markers. Among 19 markers common to both maps, the marker order on the map was maintained perfectly. CONCLUSION: Dissociation curve analysis is reliable and efficient for radiation hybrid panel scoring, and is more sensitive and robust than conventional gel-based typing methods. Several markers could be scored only by the new method, and ambiguous scores were reduced. PCR-based dissociation curve analysis decreases both time and resources needed for construction of radiation hybrid panel marker maps and represents a significant improvement over gel-based methods in any species.

The contribution of transposable elements to Bos taurus gene structure.

In an effort to identify the contribution of TEs to bovine genome evolution, the abundance, distribution and insertional orientation of TEs were examined in all bovine nuclear genes identified in sequence build 2.1 (released October 11, 2005). Exons, introns and promoter segments (3 kb upstream the transcription initiation sites) were screened with the RepeatMasker program. Most of the genes analyzed contained TE insertions, with an average of 18 insertions/gene. The majority of TE insertions identified were classified as retrotransposons and the remainder classified as DNA transposons. TEs were inserted into exons and promoter segments infrequently, while insertion into intron sequences was strikingly more abundant. The contribution of TEs to exon sequence is of great interest because TE insertions can directly influence the phenotype by altering protein sequences. We report six cases where the entire exon sequences of bovine genes are apparently derived from TEs and one of them, the insertion of Charlie into a bovine transcript similar to the zinc finger 452 gene is analyzed in detail. The great similarity of the TE-cassette sequence to the ZNF452 protein and phylogenetic relationship strongly suggests the occurrence of Charlie 10 DNA exaptation in the mammalian zinc finger 452 gene.

Comparative analysis of noncoding sequences of orthologous bovine and human gene pairs.

Genomic sequence comparison across species has enabled the elucidation of important coding and regulatory sequences encoded within DNA. Of particular interest are the noncoding regulatory sequences, which influence gene transcriptional and posttranscriptional processes. A phylogenetic footprinting strategy was employed to identify noncoding conservation patterns of 39 human and bovine orthologous genes. Seventy-three conserved noncoding sequences were identified that shared greater than 70% identity over at least 100 bp. Thirteen of these conserved sequences were also identified in the mouse genome. Evolutionary conservation of noncoding sequences across diverse species may have functional significance, and these conserved sequences may be good candidates for regulatory elements.

Comparative analysis of noncoding sequences of orthologous bovine and human gene pairs

Genomic sequence comparison across species has enabled the elucidation of important coding and regulatory sequences encoded within DNA. Of particular interest are the noncoding regulatory sequences, which influence gene transcriptional and posttranscriptional processes. A phylogenetic footprinting strategy was employed to identify noncoding conservation patterns of 39 human and bovine orthologous genes. Seventy-three conserved noncoding sequences were identified that shared greater than 70% identity over at least 100 bp. Thirteen of these conserved sequences were also identified in the mouse genome. Evolutionary conservation of noncoding sequences across diverse species may have functional significance, and these conserved sequences may be good candidates for regulatory elements.

A radiation hybrid map of bovine X chromosome (BTAX).

We present a comprehensive radiation hybrid map of the bovine X chromosome (Chr) containing 20 new markers, including both microsatellites and expressed genes. This study was conducted with a 5000-rad whole genome RH cell panel consisting of 90 hybrid cell lines. Retention frequencies of individual markers range from 7.8% for XIST to 31.1% for TGLA325. Statistical analysis with RHMAPPER placed all the loci into five linkage groups under a LOD score criterion of 6.0. These groups could be oriented relative to each other because they included multiple microsatellite loci from the consensus linkage map of the X Chr. Markers included in both this RH map and the bovine cytogenetic map were in a consistent order. The comparative bovine-human map thus generated consists of five blocks of genes, the order of which is conserved, although in the opposite direction when presented as ideograms with p and q arms. Inversions of three blocks account for the difference in gene order across the entirety of the two X Chrs.